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Four-winged flapping flyer in forward flight 

R. Godoy-Diana*’*, R Jain'’^, M. Centeno*’^, A. Weinreb^ and B. Thiria^ 


Abstract We study experimentally a four-winged flapping flyer with chord-wise 
flexible wings in a self-propelled setup. For a given physical configuration of the 
flyer (i.e. fixed distance between the forewing and hindwing pairs and fixed wing 
flexibility), we explore the kinematic parameter space constituted by the flapping 
frequency and the forewing-hindwing phase lag. Cruising speed and consumed elec¬ 
tric power measurements are performed for each point in the (/, <p) parameter space 
and allow us to discuss the problem of performance and efficiency in four-winged 
flapping flight. We show that different phase-lags are needed for the system to be 
optimised for fastest flight or lowest energy consumption. A conjecture of the un¬ 
derlying mechanism is proposed in terms of the coupled dynamics of the forewing¬ 
hindwing phase lag and the deformation kinematics of the flexible wings. 


1 Introduction 


Flapping flyers display an extremely rich variety of maneuvers because of the mul¬ 
tiple kinematic parameters that rule the unsteady production of aerodynamic forces. 
From a biological point of view, the case of four-winged flyers capable of out-of- 
phase motion between forewings and hindwings such as dragonflies is particularly 
interesting. In the words of Azuma et al. (1985): “Dragonflies can hover, fly at high 
speed and maneuver skillfully in the air in order to defend their territory, feed on live 
prey and mate in tandem formation”. Their body and wing kinematics have been 
studied extensively (Alexander, 1984, 1986; Azuma & Watanabe, 1988; Riippell, 
1989; Wakeling & Ellington, 1997a,b) and flow visualizations in tethered- and free- 
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flight configurations have demonstrated the crucial role of unsteady mechanisms 
such as the formation of leading-edge vortices in the production of lift (Thomas 
et al., 2004). Forewing-hindwing phase-lag has been shown in hovering configura¬ 
tions to be determinant for flight performance (Maybury & Lehmann, 2004): opti¬ 
mal efficiencies have been found for out-of-phase beating whereas in-phase motion 
of forewings and hindwings has been shown to produce stronger force (Wang & 
Russell, 2007; Usherwood & Lehmann, 2008). The physical mechanisms behind 
these differences in performance have nonetheless not yet been completely eluci¬ 
dated, and open questions remain in particular when considering the role of wing 
elasticity. Wing deformation is important because it can passively modify the ef¬ 
fective angle of attack of a flapping wing, determining thus its force production 
dynamics. 

In the present paper we address this problem experimentally using a four-winged 
self-propelled model mounted on a “merry-go-round”. The setup is a modified ver¬ 
sion of the one used by Thiria & Godoy-Diana (2010) and Ramananarivo et al. 
(2011), where the thrust force produced by the wings makes the flyer turn around a 
central axis. A constant cruising speed is achieved for a given wingbeat frequency 
when the thrust generated is balanced by the net aerodynamic drag on the flyer. 
These previous works have shown that passive mechanisms associated to wing flex¬ 
ibility govern the flying performance of a flapping wing flyer with chord-wise flexi¬ 
ble wings. These determine for instance that the elastic nature of the wings can lead 
not only to a substantial reduction of the consumed power, but also to an increment 
of the propulsive force. Here we introduce a new parameter using a model with two 
pairs of wings. In addition to the flapping frequency and wing flexibility, the thrust 
production is now also determined by the phase lag (p between the forewings and 
the hindwings. 


2 Experimental setup 

Figure 1 shows a sketch of the experimental setup. As in Thiria & Godoy-Diana 
(2010), the stroke plane is parallel to the shaft linking the flyer to the central bear¬ 
ing of the wheel. In addition to the four-winged instead of two-winged flyer, the 
setup also differs from Thiria & Godoy-Diana (2010) in that a counterweight has 
been added using an opposite radial shaft to balance the system. The fluctuating lift 
force is thus directed radially and absorbed by the shaft. The two wings are driven 
by a single direct-current motor with a set of gears that allows to fix the phase dif¬ 
ference between the forewings and the hindwings. All wings beat thus at the same 
frequency which was varied between 15 and 30 Hz. We have reduced the parame¬ 
ter space in the experiments reported here by fixing the physical characteristics of 
the flyer. Namely, the distance between the wings d, the stroke amplitude 9o and 
the chord-wise flexural rigidity of the wings. It should of course be noted that these 
parameters in the present tandem wing configuration, in particular the wing spacing 
d, should in general be analysed simultaneously with the forewing-hindwing phase 
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Fig. 1 Top: sketch of the experimental setup. Bottom: photos of the flapping flyer with forewing¬ 
hindwing phase lags 0 (left) and 7t (right). Each wing has a semi-circular shape so that c = 30mm 
and L = 60mm. The distance between the wings d=l mm, so that the trailing edge of the forewing 
and the leading edge of the hindwing almost touch when the wings are aligned without bending; 
the stroke amplitude max{0^„,) = max(0/,„) — do = 37°. Wings are made of 0.05mm thick Mylar 
that gives a flexural rigidity B — 3.3 x 10“^ N.m. The leading edge is thicker (1mm) and made of 
fiberglass so that it can be considered rigid and the defonnation exclusively chord-wise. 


lag (Maybury & Lehmann, 2004; Rival et al., 2011). The motion of the wings is de¬ 
scribed using the angles of the forewing and hindwing leading edges to the xz-plane, 
9fw and 9h„, respectively (see Fig. 1), as 

Sfw = 9osin(2nft) and 9h„ = 9osm{2nft — <p) , (1) 

where / is the flapping frequency and the phase lag <p is varied between 0 and iTt. 
For 0 < <p < TT the forewing is leading whereas for ;r < ^ < 2;r it is the hindwing 
that leads. The Reynolds number Re = Uc/v based on the cruising speed and the 
chord length was in the range of 1000 to 4000. 

The measured quantities are the cruising flight speed U and the consumed power 
Pj. In the study of the two winged flyer of Thiria & Godoy-Diana (2010) and Ra- 
mananarivo et al. (2011) an additional independent setup was used to measure the 
thrust force F by holding the flyer in a stationary position. The product of the force 
and cruising speed measurements was then an estimate of the aerodynamic power 
Pa = FU. One of the disadvantages of that procedure was that the two measure- 
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Fig. 2 Cruising speed and consumed power as a function of the forewing-hindwing phase lag for 
two different flapping frequencies. 


ments did not correspond to the same flight configuration: while U corresponds to 
self-propelled cruising flight, F measured at a fixed station corresponds to a ’’hov¬ 
ering” flight configuration. Here we avoid that problem by using an estimate of 
aerodynamic power obtained only from the velocity measurements, as explained in 
the next section. 


3 Results 

(i) Flying performance. The aerodynamic interactions are thus mled by (p and their 
effect can be directly measured in the performance parameters of the experiment: 
the cruising speed U and the consumed power Pj. In order to get a clear picture of 
the effect of each parameter we show first in Fig. 2 two data series corresponding to 
different flapping frequencies, plotting U and P, as a function of It can be seen 
that the phasing between the wings produces a net effect in performance, the fastest 
cruising flight velocities corresponding to a range around in-phase flapping {(p = 0), 
but the picture becomes richer when looking at the consumed power. Indeed, the 
latter has two peaks around (p « 0 and n. The previous series lie in a regime where 
increasing the flapping frequency shifts the curves to higher flying speeds and higher 
consumed power. The observed trend is readily explained using the non-dimensional 
expressions p = Pic/BcOf and u = UlAoCOf defined by Ramananarivo et al. (2011), 
where Aq is the amplitude of oscillation of the leading edge at mid span given by 
Aq = (L/2)sin0o and (Of = 2nf. The insets in Fig. 2 show the behaviour of the 
dimensionless quantities. 

The increase of speed with increasing flapping frequency is not indefinite how¬ 
ever, an effect shown in Fig. 3, where U is plotted in coloured contours in a {(p,f)- 
space for (p G [0, 7t] and scanning the full range of flapping frequencies available 
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Fig. 3 Cruising speed, consumed power and efficiency (see text) as a function of the forewing¬ 
hindwing phase lag and the flapping frequency. Only the first half o the phase-lag (p range where 
the forewing is leading was examined here. 


experimentally. For all phase lags, a clear maximum of the attained cruising speed 
occurs always around 24Hz. We will discuss in the following that this optimal fre¬ 
quency is related to the elastic properties of the wings. The second plot in Fig. 3 
shows the consumed power P,- in the same parameter space. Here the main ob¬ 
servation is that, while not surprisingly consumed power increases monotonically 
with increasing flapping frequency, the effect of the phase lag on energy expendi¬ 
ture shown previously in Fig. 2 is clearly present regardless the flapping frequency. 
Because of this effect, for different phase lags the ranges of frequency explored 
changed, giving for instance a maximum frequency for ^ = 0 of around 30Hz while 
at ^ S [;r/2,3;r/4] the frequency could reach 35Hz. We use those two measure¬ 
ments to define the following expression of efficiency, considering that the aerody¬ 
namic thrust power is proportional to (velocity times thrust force, the latter being 


1 = 


\pSU^ 

Pi 


( 2 ) 


where S is the effective wing surface. Other definitions of efficiency using purely 
dynamical parameters (e.g. Kang et ah, 2011) should give equivalent results to the 
expression 2 chosen here in terms of the measured consumed power P, . It can be seen 
that the optimum in terms of efficiency is shifted toward larger phase lags (around 
(p « 7t/2) than the optimum in terms of maximum cruising speed. 


(ii) Wing kinematics. In addition to the performance measurements, the wing mo¬ 
tion was recorded using high-speed video in a fixed frame. That is, with the mechan¬ 
ical insect not mounted on the merry-go-round but on a hxed base. Figures 4 and 5 
show typical time series for (p = 0 and 1.37r, respectively. A dark screen was used to 
mask the right side wings and have a clearer view of the kinematics. The main fea- 
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Fig. 5 Sequence showing the kinematics of the flapping wings at / = 23Hz and (p = l.3n 
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Fig. 6 Tracking the wing motion. Example of in-phase flapping at / = 23Hz. Black circles are 
the points tracked on each wing for the leading and trailing edges (FWLE and EWTE, correspond 
to the forewing, while FIWLE and FIWTE correspond to the hindwing). The open circles are the 
positions of the win tips tracked. 


ture that can be seen in these time series of snapshots is the large deformation of the 
wings during the flapping cycle. Indeed, the wings bend over a length scale that is of 
the same order of magnitude than the chord length. In addition to the kinematics of 
the compliant wings, which can be followed for each wing independently, figures 4 
and 5 hint on the complex interaction that arises from the combination of the bend¬ 
ing dynamics and the imposed forewing-hindwing phase lag. We will analyse this 
interaction in the following by tracking the motion of different points on each wing. 

Figure 6 presents the tracking in time of the (x,y) positions of the leading and 
trailing edges of each wing at the point of maximum chord length (black circles, 
FWLE and FWTE denoting forewing leading and trailing edge, respectively, while 
FIWLE and HWTE correspond to the equivalent points on the hindwing). Addi¬ 
tionally, the tip of each wing is also tracked (empty circles). The latter is useful to 
minimize the measurement error in the leading edge amplitude of motion since it 
has a larger swept amplitude and the relative error is thus smaller. The time-series 
of the positions of these points are shown in figure 7 (a) and (b) for two different 
values of (p. The main observation here is the trailing-edge-leading-edge phase lag 
7 for each wing, which has been reported in Ramananarivo et al. (2011) to be a 
crucial element of the propulsive performance of flexible flapping wings. The mea¬ 
sured value of 7 for several forewing-hindwing phase lags is shown in figure 7 (c). 
A slight decrease is observed for 7 when (p decreases from 2n to L3;r, i.e. when the 
systems goes from in-phase flapping to a configuration with the hindwings leading 
by approximately a third of a period, but the effect is very weak. 
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Fig. 7 Time series of the j'-position of the leading and trailing edges of both wings for / = 23Hz 
and (a) (p = 0 and (b) l.37t. The leading-edge-trailing-edge phase lag 7 is shown schematically. 
In the 9 = 0 case, 71 = 75 = 7 . The y-axis is rendered dimensionless using A, the peak-to-peak 
amplitude swept by the leading edge at mid-chord, (c) Trailing-edge-leading-edge phase lag 7 for 
different forewing-hindwing phase lags. The two different markers correspond to the forewing 71 
(A) and hindwing (v). 



4 Discussion 

Before discussing the effect of the forewing-hindwing phase lag we start with a 
comment on the role of the flapping frequency. The existence of an optimal flapping 
frequency for the present setup (as is evident in Fig. 3) is related to the effect of 
wing compliance on the propulsive performance of the flapping wings, a question 
that has been widely discussed in the literature (Shyy et al., 2010; Spagnolie et al., 
2010; Masoud & Alexeev, 2010; Thiria & Godoy-Diana, 2010; Ramananarivo et al., 
2011; Kang et al., 2011). For flexible wings flapping in air, where the mass ratio of 
the wing with respect to the surrounding fluid is high, the main bending motor is the 
wing inertia (Daniel & Combes, 2002; Thiria & Godoy-Diana, 2010). Increasing the 
flapping frequency leads to an increased deformation of the wing, which is useful 
in terms of propulsive performance up to a certain point where the effective lifting 
surface is diminished and thrust production drops (Ramananarivo et al., 2011). An 
interesting point is that the optimal frequency that can be estimated by inspec¬ 
tion of Fig. 3 is different depending on whether one considers maximum cruising 
speed (fopt 24 Hz) or maximum efficiency {fopt 21 Hz). This can be explained 
by the fact that the consumed power, which enters the quotient defining efficiency 
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(Eq. 2) in the denominator, increases monotonically with an increasing frequency. 
This will thus shift the location of the peak in the efficiency that corresponds to the 
aerodynamic power peak. 

In order to focus on the effect of the forewing-hindwing phase lag, a fixed fre¬ 
quency can be chosen. In Fig. 2, for each of the two frequencies represented, the 
same trends are observed in the cruising speed and the consumed power and, as 
shown in the insets, the frequency dependence is well explained using the “elasto- 
inertial” dimensionless variables u and p defined above. Now, with respect to the 
phase lag (p, two main points appear clearly: on the one hand, the fastest flying per¬ 
formance is found around in-phase flapping (i.e (p — 0), while around anti-phase 
flapping {(p = n) this cruising speed is the lowest. We may note that the curve is 
not symmetric with respect to ^ = 0, since it diminishes only slightly (one could 
consider a plateau) between (p — 0 and njl, whereas it drops more rapidly when 
the hindwing leads (i.e. when going from (p = 0 = 27t towards (p = 3;r/2). The con¬ 
sumed power curve on the contrary has two maxima (around (p —0 and n) and two 
minima around ((p — ;r/2 and 3;r/2). In-phase and anti-phase flapping being the 
most expensive can be explained by a simple inertial argument since during these 
configurations the motor has to accommodate the acceleration/deceleration of both 
pairs of wings at the same time, contrary to the intermediate phase lags (p ^ n/2 
and 3n/2 where the deceleration of one pair of wings occurs while the other pair is 
accelerating, hence redistributing the load on the motor. 

The shape of the cruising speed and consumed power curves determines the lo¬ 
cation of the maximum observed in the efficiency contours in Fig. 3 being around 
(p — nl2. A secondary maximum can be expected around ^ = 37r/2 (a zone of the 
parameter space that was not fully explored in the present experiments). We can 
now compare the optimum phase-lags that lead to peaks of cruising speed and of ef¬ 
ficiency and comment on their physical origin: while the maximum cruising speed is 
observed for phase-lags between zero and ;r/4, the optimum phase-lags in terms of 
efficiency are around 7r/2. The former are ruled solely by the aerodynamics, where 
the performance of different kinematic patterns will have to be analysed considering 
for instance the interaction between hindwing and the vortex structures shed by the 
forewing. A wake capture process of this sort has been proposed by Kolomenskiy 
et al. (2013) as a possible explanation for the large propulsive force found at Q.lSn 
in their 2D numerical simulations.* Concerning the optimum phase lag in terms of 
efficiency requires on the other hand considering the power consumption, which 
is not only correlated to aerodynamics, but has a large contribution determined by 
solid inertia, as we have mentioned in reference to the power plot in Fig. 2b. 

The effect of the modulation of (p in terms of aerodynamics, as mentioned pre¬ 
viously, will be intrinsically related to the roles of the distance d between the two 
wing pairs and of the deformation kinematics. In this paper we have fixed d and 
considered a single flexibility in order to explore the possible roles of the wing de¬ 
formation in the propulsive performance. Considering our previous studies on the 
effect of flexibility with a two wing flyer (Thiria & Godoy-Diana, 2010; Ramana- 

* Note that in Kolomenskiy et al. (2013) the phase lag is defined with a negative sign with respect 
to (p used here so that O.lSll here corresponds to their 1.25;r. 



10 


Godoy-Diana et al. 


narivo et al., 2011 ), where the trailing-edge-leading-edge phase lag 7 was shown 
to be a crucial parameter to determine performance, we attempted here to monitor 
modifications in 7 as a function of (p (see the analysis of wing kinematics in Figs. 

6 and 7). We picked a range of (p (decreasing from 2n to 1.3;r) where a clear de¬ 
terioration of the flying performance is observed in figure 2 as the hindwing starts 
leading the forewing. The observations on the trailing-edge-leading-edge phase lag 

7 are however not conclusive. A slight decrease of 7 is indeed observed while (p 
goes from 2n to 1.3;r, which at these frequencies could explain a diminishing thrust 
performance, but the effect is too weak and the present results do not permit to give 
a thorough and quantitative confirmation of the effect if it exists. 


5 Conclusions 

We have shown that a four-winged flapping flyer in a cruising regime does present 
different optimal forewing-hindwing phase lags depending on whether one would 
want to tune the system for maximum cruising speed or minimum energy expendi¬ 
ture. These results are in accordance with previous studies in hovering configura¬ 
tions by Wang & Russell (2007); Usherwood & Lehmann (2008), hinting that the 
mechanisms described here should be robust elements to consider in any aerody¬ 
namic model: one the one hand wing inertia is a major player in the power expen¬ 
diture, while on the other hand the thrust production is ruled by the aerodynamics 
around the flexible wings. 

A full flow field reconstruction around the wings is certainly desirable to compare 
different points in the parameter space with different performances and clearly iden¬ 
tify the aerodynamic mechanisms at play. Experimentally this can be challenging, 
and numerical simulations such as a 3D version of the fluid-structure simulations of 
Kolomenskiy et al. (2013) can be an interesting option to define robust aerodynamic 
models. 

Other issue that remains not fully understood is the effect of the forewing¬ 
hindwing separation d as a parameter independent of the forewing-hindwing phase 
lag which will have a different role depending on whether the system is hovering or 
cruising. This point will be of particular importance when going beyond the “steady” 
regimes of hovering or cruising and into the study of transient regimes. These bring 
indeed a vast set of open questions related to the unsteadiness —like for instance the 
multi-body dynamics of manoeuvring— and where the wings are to be analysed as 
part of a full system accelerating, performing sharp turns (e.g. Bergou et al., 2010) 
or taking off (e.g. Bimbard et al., 2013). 
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